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We estimate the 7 = scalar meson <r//o(600) parameters from 7T7t and 77 scattering data below 700 MeV using an improved 
analytic K-matrix model. A fit of the hadronic data gives a complex pole mass M a = 422 — i 290 MeV, while simultaneous 
best fits of the 77 — > 7r + 7r _ , 7r°7r° data give a direct width of (0.13 ± 0.05) keV, a rescattering component of (2.7 ± 0.4) keV 
and a total (direct +rescattering) width of (3.9 ± 0.6) keV. "Running" these results to the physical real axis, the small "direct" 
77 and the large hadronic widths at the "on-shell" mass are compatible with QCD spectral sum rules (QSSR) and some low- 
energy theorems (LET) expectations for an unmixed lowest mass glueball/gluonium as of a mass around 1 GeV and a large 
OZI-violation decay into nn . 



1. Introduction 

Understanding the nature of scalar mesons in terms of 
quark and gluon constituents is a long standing puzzle in 
QCD lj. One might expect that the decay rate of these 
mesons into two photons could provide an important 
information about their intrinsic composite structure. 
The problem here is that some states are very broad {a 
and k mesons), others are close to an inelastic thresh- 
old (/ (980), a (980)), which makes their interpretation 
more difficult. Besides the interpretation within a qq 
model [1I2I3I4I5I6I7] or unitarized quark model |8l9j , also 
the possibility of tetraquark states |10lllll2ll3ll4j (and 
some other related scenarios: meson-meson molecules 
|15ll6j . meson exchange [17]) is considered. In addi- 
tion, a gluonic meson is expected in the scalar sector, 
according to lattice QCD [18|19j . QCD spectral sum 
rules (QSSR) [20l2ll22l23l24T4T7] a la SVZ g Ug], and 
some low-energy theorems (LET) |24I27I28| . Such a 
state could mix with the other qq mesons I3|4)6| 29(30 ] . 
Among the light particles, the /o(600)/cr meson could 
be such a gluonic resonance. It can manifest itself in 
some effective linear sigma models [31 32 33 or con- 
tribute to the low-energy constants at 0(p 4 ) of the 
QCD effective chiral Lagrangian [34]. Recent analy- 
ses of the 77 — > 7T7T processes have extracted the width 
of / (600)/cr -> 77 in the range: (4.1 ± 0.3) keV [35], 
(3.5 ±0.7) or (2.4 ±0.5) keV [36 to (1.8 ±0.4) keV [37], 
while the one from nucleon electromagnetic polarizabili- 
ties has given (1.2 ±0.4) keV [38] . Some of these results 
have been interpreted in [35] as disfavouring a gluonic 
nature which is expected to have a small coupling to 77 
[241314114116113] . In the following, we shall reconsider 
the analysis of the 77 — > tttt process in the low energy 
region below 700 MeV, where we conclude that it is 
dominated by the coupling of the photons to charged 
pions and their rescattering, which therefore can hide 
any direct coupling of the photons to the scalar mesons. 
Some first results from this study have been presented 
in [Ml- 

* Email: gerard.mennessier@lpta.univ-montp2.fr 
t Email: snarison@yahoo.fr 

* Email: wwo@mppmu.mpg.de 



2. Results suggesting a light ++ glueball 

The existence of glueballs/gluonia is a characteristic pre- 
diction of QCD and some scenarios have been developed 
already back in 1975 [3D]. Today, there is agreement 
that such states exist in QCD and the lightest state has 
quantum numbers J PC = ++ . 

Lattice QCD - Calculations in the simplified world 
without quark pair creation (quenched approximation) 
find the lightest state at a mass around 1600 MeV [18] . 
These findings lead to the construction of models where 
the lightest glueball/gluonium mixes with other mesons 
in the nearby mass range at around (1300-1800) MeV 
(see, for example, [29]). However, recent results beyond 
this quenched approximation [19] suggest that the light- 
est state with a large gluonic component is expected in 
the region around 1 GeV, and therefore, a scheme based 
on the mixing of meson states with all masses higher 
than 1300 MeV could be insufficient to represent the glu- 
onic degrees of freedom in the meson spectrum. Further 
studies concerning the dependence on lattice spacings 
and the quark mass appear important. 

QSSR and LET - These approaches have given quan- 
titative estimates of the mass of glueballs/gluonia [2"2l 
12314] and of some essential features of its branching ra- 
tios [241417120131] , As we shall see in section[S] these ap- 
proaches require the existence of a gluonium erg below 1 
GeV having a large 7T7t but small 77 widths, in addition 
to its corresponding radial excitation a' B and a narrow 
gluonium G(l. 5-1.6) coupled weakly to tttt [221241417] . 

Phenomenological studies - The identification of the 
scalar glueball from experiment requires a full under- 
standing of the scalar meson spectrum. After group- 
ing states into flavour multiplets the left over states are 
candidates for glueballs. Several schemes exist, moti- 
vated by the quenched lattice result, where the extra 
gluonic state is assumed to mix into the three isoscalars 
/o(1370), / (1500) and /n(1710) [29]. In an alternative 
approach [6] (see also ref. [41]) the lightest qq mul- 
tiplet is formed from / (980), / (1500), 7^(1430) and 
ao(1450), and the glueball is identified as the broad ob- 
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ject at smaller mass represented by Jo (600; it dominates 
7T7T scattering near 1 GeV but extends from threshold up 
to ~ 1800 MeV. The appearance of this broad object in 
most gluon rich processes was considered in support of 
this hypothesis. In this analysis of the spectrum, results 
from elastic and inelastic tttt scattering as well as from 
D, B and J/tp decays have been considered |6I42I43| . 

3. Constraints on 77 from tttt processes 
Given a set of coupled multi-channel strong processes, 
related electromagnetic or weak processes are largely 
predetermined. First, there are the multi-channel ex- 
tended unitarity relations for the amplitudes with mixed 
strong and weak couplings (corresponding to the Watson 
theorem for the single channel); secondly, there are dis- 
persion relations which these amplitudes have to fulfill 
(though not proved for composite particles). The gen- 
eral consequences of these constraints have been inves- 
tigated by Omnes and Muskhelishvili [44 for the single- 
and multi-channel cases. The amplitudes for the weakly 
coupled processes are largely determined by strong am- 
plitudes, but there are polynomial ambiguities which re- 
sult from subtraction terms of the dispersion relations 
depending on the asymptotic high energy behaviour. 

• The analytic K-matrix model 

This general formalism has been applied by Mennessier 
[45] (see also [46]) to the calculation of the electromag- 
netic processes 77 — ► tttt, KK, given the strong pro- 
cesses tttt — > 7T7T, KK. The strong processes are repre- 
sented by a K matrix model representing the amplitudes 
by a set of resonance poles. In that case, the dispersion 
relations in the multi-channel case can be solved explic- 
itly, which is not possible otherwise. This model can 
be reproduced by a set of Feynman diagrams, including 
resonance (bare) couplings to tttt and KK and 4-point 
tttt and KK interaction vertices. A subclass of bubble 
pion loop diagrams including resonance poles in the s- 
channel are resummed (unitarized Born). The model 
also includes contributions from the exchange of vec- 
tor mesons in the t-channel which become important 
at the higher energies above 700 MeV. The 77 scatter- 
ing amplitudes also fulfill the constraints at s — re- 
quired by the Born approximation. The radiative width 
of the resonances cannot be predicted due to the poly- 
nomial ambiguity, which can be taken into account by 
introducing as free parameters the "direct couplings" 
of the resonances to 77 in an effective interaction ver- 
tex. In the present analysis we have extended the model 
by the introduction of a shape function which takes ex- 
plicitly into account left-handed cut singularities for the 
strong interaction amplitude. This allows a more flexi- 
ble parametrisation of the tttt data at low energies and 
improves the high energy behaviour. Next we discuss 
the features of the model at low and at high energies. 

• Low energy limit: rescattering 

A striking feature of the low energy 77 — > tttt scattering 
is the dominance of the charged over the neutral tttt cross 
section by an order of magnitude which can be explained 
by the contribution of the one-pion-exchange Born term 



in 77 — > tt + tt~ . In the process 77 — > tt°tt , the photons 
cannot couple "directly" to 7r°7r° but through interme- 
diate charged pions and subsequent rescattering with 
charge exchange. This feature is realized in the analytic 
model [45 . A similar approach has been applied in the 
model of Ref. [47] which described data on 77 — * tttt up 
to 900 MeV by rescattering. A modern example for the 
importance of rescattering is Chiral perturbation theory. 
To one-loop accuracy, the calculation based on one pion 
exchange and tttt rescattering relates the 77 — > tt°tt 
to the tt + tt~ —* tt°tt cross-section using m n , /„• and 
e [48 . The 2-loop corrections obtained in [49] includes 
contributions from other exchanges (spin 1). It amounts 
to about 30% of the lowest order results and provides a 
reasonable description of the data up to ~ 700 MeV. 

• High energy limit: direct production 
Whereas at low energy the photon interacts by its cou- 
pling to the charged pions, the situation becomes differ- 
ent at high energies where the photon can resolve the 
constituents of the hadrons. An example is the produc- 
tion of /2(1270) in 77 — > tttt. The analytic model [45] 
predicts a decreasing cross section for this process with 
pion exchange reaching < 10% of the observed cross sec- 
tion at the peak of the f%. After inclusion of vector ex- 
change, this contribution doubles. Because of form fac- 
tor effects involved in the photon coupling to hadrons, 
these predictions from point like hadrons are to be taken 
rather as upper limits. This requires the need of "direct 
coupling" of f% — * 77. Indeed, it is well known that 
the radiative decays of the tensor mesons f^, o-i are 
well described by a model with direct coupling to the 
qq constituents according to the SU(3) structure of the 
nonet with nearly ideal mixing (see, e.g. Ref. 50\). An- 
other example of a direct process is the production of p 
mesons in the process 7p — » tt + tt~p by direct VMD cou- 
pling, which is found much larger than the rescattering 
contribution from the "Drell-Soding" background pro- 
cess [51]. We therefore interpret the direct terms in the 
77 processes as representing the amplitudes with cou- 
pling of the photon to the partonic constituents of the 
resonances (such as qq, 4q, gg,. . .). 

4. Direct and rescattering couplings of the a 

In the application of the analytic model [45], we first 
obtain a suitable parametrisation of the 7T7t scattering 
data, and then determine the direct coupling by com- 
paring the model with the 77 results. In the present 
Letter, we restrict ourselves to the low mass region be- 
low 700 MeV where we neglect vector and axial- vector 
exchanges 0, inelastic channels and D- waves, further- 
more, we assume a pointlike pion-photon coupling. 

• Amplitude for elastic tttt scattering 

In the K-matrix fits to the tttt elastic scattering data for 
the range of energy [(0.5-1.8) GeV] , a pole is found in the 
isoscalar S-wave amplitude with a large imaginary 
part which corresponds to a state |52|53|54] : 
M ps T ps 1 GeV, (1) 

4 They only start to be relevant above 700 MeV |i5) . 
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though this value may depend on the treatment of the 
other resonances 13016141] . The mass value in Eq. ([1]) 
is close to the one where the observed S-wave phase 
shift goes through 90° as in a simple Breit Wigner form 
without background. This broad object has been called 
/o(1000) [Si], for a while / (400 - 1200) by the PDG 
and now /o(600)/er. In general, the complex resonance 
self-energy £ = M — (i/2)T is energy dependent and 
determines the poles of the S-matrix corresponding to 
mass and width of the physical particles; they do not 
usually coincide with the Breit Wigner masses. In recent 
determinations, where analyticity and unitarity proper- 
ties are used to continue the amplitude into the deep 

imaginary region, one obtains for the compl ex a poleH: 
441 - i 272 MeV [SS] or 489 - i 264 MeV [SSj . (2) 

We shall see that the a pole can be referred to the 
same broad object defined above. For applying the an- 
alytic model [45j . we introduce a shape function fo(s) 
which multiplies the cnnr coupling. For simplicity, we 
do not include the 4-point coupling term. The real ana- 
lytic function fo(s) is regular for s > and has a left cut 
for s < 0. For our low energy approach, a convenient 
approximation, which allows for a zero at s = sa and a 
pole at ctd > simulating the left hand cut, is: 

S - SAO 



Ms) = 

S + (JDO 

The unitary tttt amplitude is then written as: 

T (o) M _ Gf (s) _Gf (s) 
1 o \ s ) 



(3) 



- Gf(s) V(s) 
■ lS i' } sintff } /P with p( s ) = (1 - 4m2/ s )i/2 ; q 



(4) 



SR 



<j£ B is the bare coupling squared and: 

lmV = lm(-Gfo) = -(6p)G f , (5) 
with: (9p)(s) — below and (9p)(s) — p(s) above 
threshold s = Am\. The amplitude near the pole sq 
where T>(sq) — and V{s) « V{so){s — sq) is: 



,(0) 



9-1 = 



G/o(s ) 



(6) 



so-s -W(so) 
The real part of T> is obtained from a dispersion relation 
with subtraction at s = and one obtains: 

■ ■ (7) 



ho(s) -ho(0) 



Ms) = - 

TT 

where: h (s) = f (s)L s i (s)-(ctno/ (s+o-do))L s i (-ctdo), 

a no is the residue of /o(s) at — <tdo and: L s i(s) — 

[ (s — 4m^) /m^]Li(s, ml) with L ± from [45] . 

• Amplitudes for 77 — > tttt scattering 

Starting from the S wave amplitude in Eq. ^ we derive 



the amplitude T 7 
isospin / = as: 



(i) 



for the electromagnetic process for 



(8) 



Here the contribution from the Born term of 77 — > 
ir + n~ is given by f$ = 2L\ as defined in [45], a real 
analytic function in the s plane with left cut s < 0. The 




function /q represents tttt rcscattcring; it is regular for 
s < 4m^ but has a right cut for s > 4mJ with: 

Imf*(s + ie) = (9pf f*)(s) , (9) 
which vanishes at s = 0. With this definition the Wat- 
son theorem is fulfilled, i.e. the phase of is the same 
as the one of the elastic amplitude I? -1 in Eq. ((4]). The 
real part is derived from a dispersion relation with sub- 
traction at s = for to satisfy the Thomson limit, 
and has a representation similar to the one in Eq. (O, 
but by replacing (fo,ho) by (/(f )■ The function 
is defined as ho below Eq. ([7]) but with L\ replaced by 



-Lf everywhere. It vanishes at 



-a D0 and ffr (s) is 



regular at this point. Finally, the polynomial V reflects 
the ambiguity from the dispersion relations and is set 
here to V = sF 7 ^/2. It represents the direct coupling of 
the resonance to 77. The residues at the pole so of the 
rescattering and direct contributions to in Eq. ((8]), 
respectively, are obtained as: 

B ' s o). _dir _„*oF 1 V2 

9ir — 



o rcsc o 



2 GM 
3 a -V'(s ) 



4 h 



from which one can deduce the branching ratio: 

2 

Fl . 



r 

*- tr— >77 
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2 2a 2 


so 


■Tcr — >1T7T 


" \P(S0)\ 




" Wo)\ 


GMso) 



(10) 



(11) 



Similarly, we parametrize the 1 = 2 S-wave amplitude 
(2) 

Tq ; by introducing the shape function / 2 : 



rp(2) _ 



A/a («) 



l-A/ 2 (s)' 
and obtain: 



T (2) 
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where /f = / f 



Ms) 



SA2 



(S + (T D l)(s + <JD2) 




-A/ 2 ( S )J ' 

: Im/ 2 (,s) = (6p)f 2 (s), Im/. 



, (12) 

(13) 
! ( S ) = 



[0p f2 / 2 B )(s)- These amplitudes are again both sub- 



12 J2 

tracted at s 
analogy: 

Ms) = - 



= as in case of / = and one finds in 



h 2 (s)-h 2 (0) 



(14) 



where: h 2 (s) = f 2 (s)L s i(s)-(a N i/(s + cr m )) L s i(-a m ) 
— (o~N2/(s + o- D 2))L sl (~a D2 ); (r Nu a N2 are the residues 
°f f2( s ) at — Uui, — (7d2 an d f 2 (s) is defined as f 2 (s) 
in Eq. (JT5J) but with L sl replaced by — L 2 sl , The cross 
sections for the tttt and 77 scattering processes are ob- 
tained from the previous expressions of the amplitudes 
as described in the appendices of [45], where the cos 9 
integration should be done from to 1 for the neutral 
and from -1 to 1 for the charged cases H. 
• Results of the analysis 

The analytic model in its original parametrization [45] 
without direct term has given already quite a decent de- 
scription of the data available now for both charge states 
7r + 7T~ and 7r°7r° within ~ 20% as discussed previously 
[39] . This already indicates that the direct term is much 
smaller than the rescattering term. In the present anal- 
ysis, we repeat the calculation but take the new version 
of the model [45] by introducing the shape junction fo(s) 



5 A similar result [see Eq. H17l ll is found in the analytic model |45| . 



6 The range of integration may have been misused in 'f>7'_ . 
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defined in Eq. ([3]). We first determine the parameters 
of the model for tttt scattering below 700 MeV from the 
best approximation of our formula to the phase shifts 
5^ and 5^ obtained in the Roy equation analysis in 
[55] , which takes into account the high-energy behaviour 

2. 



suggested in [5j>] . For / 

sao = 0.016, 
s R = 0.823, 

and for 1 = 2: 
s A 2 = 0.039, 
ctdi = 0.693, 



= 0, we obtain in units of GeV 
G ~ 1.184 , 
arm = 0.501 , 



A ~ -4.136 , 
o~ D2 = 4.384 , 



(15) 



(16) 



(17) 



from which, we find the pole mass and residue: 
M a ~ 422 - i 290 MeV ; 57r ~ 0.06 - i 0.50 GeV, 
which is close to the mass M„ in Eq. and in the 
previous paper [35]. Given the tttt amplitude we can 
predict the cross sections for 77 — * tttt where the only 
free parameter F 1 is related to the strength of the di- 
rect coupling a — > 77. The fit of the model involving 
both the tttt rescattering and direct meson couplings to 
the data from Crystal Ball (7r°7r°) [58J and MARK-II 
(tt + tt~) [59] collaborations is shown in Figs. [T] and [2] for 
the cases with and without direct contribution. A si- 
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Fi gure 1. Fit of the 7r^7r^ cross-section in nanobarn (nb) versus 
^fs using unitarized Born amplitude: F 1 = (dot-dashed); F 1 = 
-0.09: 1=0 (large dashed), 1=0+2 (continuous); F 7 = -0.07: 
1=0+2 (small dashed). The data are from Crystal Ball I58| for 
|cos#|<0.8; 

multaneous best fit of the data in the region below 700 
MeV, is obtained for : F 1 ~ —0.070, at the minimum 
X 2 — 25.6 (19 data points) 0. However, a fit of the neu- 
tral channel alone gives a much better fit: F 7 ~ —0.090, 
for a x 2 — 3-38 (8 data points). In the charged channel, 
this is due to the deviations of the fit at high-mass; the 
large systematic errors and the absence of data points 
in the region 0.40 to 0.55 GeV do not permit a good 
understanding of this channel. One can notice that the 
1 = 2 S'-wave amplitude which reproduces with high- 
accuracy the CGL points [55) improves the agreement 
with the data in the neutral case (Fig. [1]) but gives 



7 In the earlier analysis |45| . a value F~ of 0(1) has been found, 
which corresponds to the DM2 data in 60 which are about twice 
the one of MARKII '59 , i.e twice the Born contribution. 
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Figure 2. The same as in Fig. [T] but for 7r + 7r~ using unita- 
rized Born amplitude: _F 7 = (dot-dashed); F 7 = —0.07 and 
1=0+2 (small dashed). The continuous line corresponds to the 
non-unitarized Born amplitude with F 1 = 0. The data are from 
MARKII [59] for | cos0| < 0.6. 

negligible contribution for the charged pion. For defi- 
niteness, we consider as a final estimate the mean value 
of the previous numbers: 



F 7 ~ -(0.080 ±0.014) , 



(18) 



where we consider that the uncertainty in the fitting pro- 
cedure 0.010 is expressed by this range. We have added 
linearly the systematics of the method to be about 5% , 
which we have estimated from the deviation of the mass 
and width determinations in Eq. (|17p from the ones 
[55] in Eq. j2]), as we have used their parametrization 
for fixing the parameters of the model. The correspond- 
ing values of the residues in units of a x 10~ 3 GeV are: 

(7±l) + i (31 ±4) , 
gpf~ ~ (90 + 5) +i (116 ±6) . (19) 

Using Eqs. ([Tl ] ITT l [15] ) and r ff ^ X7r = 21mA/,,, we can 
deduce the "partial" 77 widths at the complex pole : 
(0.13 ±0.05) keV , 

(2.7 ± 0.4) keV , (20) 



pair 

a — >77 
presc 

a — >77 



and the total 77 width (direct 

ptot 



rescattering): 



iT . ^ (3.9 ±0.6) keV , (21) 
Improvements of these estimates need more precise data 
below 700 MeV, and an extension of the analysis to 
higher energies. 

• Comparisons with existing results 

Direct 77 width - Our result for the direct coupling is 
comparable with our previous value [39] obtained with- 
out using the shape function fo(s). In a parallel analysis 
[61], a K-matrix model analogous to Ref. [45] has been 
applied in a larger energy range of energy, including the 
recent BELLE data [55] and still assuming elementary 
pion and kaon exchange; a smaller direct coupling was 
obtained, but without error analysis. 

Total 77 width - Our result is compatible with the 
range of values (1.2 ~ 3.2) keV obtained in [39] using a 
Breit-Wigner parametrization of the 77 — > tt°tt data. 
It is also in the range of values based on dispersion rela- 
tions (see introduction), though a direct comparison is 
difficult to do as these analyses are extended to higher 
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energy region, and we have an additional determination 
of the direct term. In our analysis below 700 MeV, we 
have also neglected the oj and axial-vector exchange in 
agreement with |45j . Other calculations find contribu- 
tions of 2.5% 05] or 16% [37] at 400 MeV. One can 
notice that the effect of the direct term in our Fig. [TJ 
distorts the shape of the unitarized Born curve, and then 
permits a good fit of all points below 700 MeV, which 
is not the case of the one in [37]. Then, if our fit was 
restricted to lower masses to suppress these other ex- 
changes no different results are to be expected. 
• "On-shell" a mass and widths 
For a more appropriate comparison with QSSR predic- 
tions , one can introduce the mass and widths of the 
"visible meson" on the real axis. This can be either the 
Breit-Wigner mass and width in Eq. ([1]) or/and the "on 
shell" mass and widths (see, e.g. [63]). 

Mass oj the a - A Breit-Wigner parametrization of 
the data, leads to value given in Eq. H}: M£ w ps ps 
1 GeV, which can be compared with the "on-shell mass" 
M° s in the model 05] where the amplitude is purely 
imaginary at the phase 90°: 

ReX>((M° s ) 2 ) = => M° s ps 0.92 GeV . (22) 

Hadronic width - In the same way as for the mass, 
one can define an "on-shell width" [see Eqs. ([5]) and ©] 
evaluated at s = (M° s ) 2 : 

M ™ T ° - =► r ™ * L02 GeV > ( 23 ) 

to be compared with the Breit-Wigner width in Eq. (I]) . 

Direct 77 width - As there are no separated S wave 
cross sections we use here only the estimates based 
on the analytic model [45] where we proceed as in 
the previous section but we evaluate the parameters at 
s = {M° s ) 2 . Then, we obtain from Eq. {JTTJ : 



(1.0 ±0.4) keV. 



(24) 



pos.dir 

a — >77 

As the model is extrapolated here towards energies be- 
yond its validity, we expect that the previous results are 
a crude approximation and will only serve as a guideline. 

5. QSSR/LET results and QCD tests of the a 

The QSSR/LET results obtained in the physical region 
can be better compared with experiment by using the 
result for a Breit-Wigner parametrization of the data in 
Eq. (fT]) or by using the results for the on-shell mass de- 
rived previously but not by directly comparing with the 
parameters of the complex pole. The direct 77 coupling, 
which can reveal the photon coupling to the intrinsic 
quark [qq, Aq, . . .) or/and gluon ( gg, . . .) constituent in- 
ternal structure of the resonance, can be also related to 
the QSSR and/or LET evaluations of its width through 
quark or gluon loops like is the case of the quark trian- 
gle for the pion and/or the /2(1270). As the comparison 
with QSSR/LET is an important step for interpreting 
our fitted values, it is essential to review shortly the 
main steps in the derivations of the results. 
• Gluonia masses from QSSR 

Masses of the bare unmixed scalar gluonium can be de- 
termined from the two Laplace unsubtracted (USR) and 



subtracted (SSR) sum rules: 
1 



£-i(t) = -V>(0) + 



dte- tT lTmjj(t) , 
dt 



1 



t 



e- fT Iim/>(<) , 



(25) 



of the two-point correlator ip(q 2 ) associated to the trace 
of the QCD energy-momentum tensor current: 
1 



-B{a s )G%G^ + [l+ 7 mK)] E ™«^VV (26) 



u,d,s 



t is the sum rule variable; f3(a s ) = 3i(a s /ir) + ... and 
7 m (a s ) = 71 (as/71") + ... are the QCD (3- function and 
quark mass anomalous dimension: — (3i = (1/2)(11 — 
2n//3), 71 = 2 for SU(3) C x SU(n f ). The subtraction 
constant V'(O) = — 16(/3i/7r)(a s G 2 ) expressed in terms of 
the gluon condensate [20] (a s G 2 ) = (0.07 ±0.01) GeV 4 
|64l26j affects strongly the USR analysis which has 
lead to apparent discrepancies in the previous litera- 
ture when a single resonance is introduced into the spec- 
tral function [23]. The SSR being sensitive to the high- 
energy region (t ~ 0.3 GeV -2 ) predicts 0]: 

Mq ~ (1.5 - 1.6) GeV , (27) 

comparable with the quenched lattice value [18] . while 
the USR stabilizes at lower energy (r ~ 0.8 GeV -2 ) and 
predicts a low-mass gluonium [22j : 

M aB ~ (0.95 - 1.10) GeV , (28) 

comparable with the unquenched lattice value |19] . Fur- 
thermore, the consistency of the USR and SSR can be 
achieved by a two-resonances (G and a) + "QCD con- 
tinuum" parametrization of the spectral function [24pl . 

• Gluonia couplings to Goldstone boson pairs 

These couplings can be obtained from the vertex func- 
tion: 

U[ <Z 2 = ( <Zl - (Z2 ) 2 ] = (7r 1 |^|7r 2 ) , (29) 
obeying a once subtracted dispersion relation: 

V{q 2 ) = V(0) + q 2 / —-TmV(t) . (30) 

J 4m 2 t t- q z -ieir 

with the condition: V(0) — Qirnfy — > in the chiral 
limit. Using also the fact that V'(0) = 1, one can then 
derive the two sum rules: 
1 
4 



E 

S = <JB,- 



V2/ s =o, - e 



S — <7 S • 



V2fs 
Ml 



1,(31) 



where fs is the decay constant analogue to The 1st 
sum rule requires the existence of at least two resonances 
coupled strongly to irir. Considering the ob and a' B but 
neglecting the small G-coupling to tttt as indicated by 
GAMS REE one predicts in the chiral limit [2414) : 

(4 - 5) GeV , (32) 



\9. 



\9c B K+K-\ 



8 In [4] the QCD continuum has also been modelized by a a' g (ra- 
dial excitation of the <tb), which enables to fix the decay constant 
f ag:tT ' once the as, <?' B masses are introduced as input. 
9 The G(l .6) is expected to couple mainly to the U(1)a channels 
r/r/ and through mixing to r?r/ , r\r\ 24 4 . 



6 



a universal coupling, which will imply a large width Lj: 

12 / a^2 \ V 2 



1.9a 



16ttM 



1 



4m; 



0.7GeV.(33) 



This large width into irir is a typical OZI-violation due 
to non-perturbative effects expected to be important in 
the region below 1 GeV, where perturbative arguments 
valid in the region of the G(1.5-1.6) cannot be applied. 
This result can be tested using lattice calculations with 
dynamical fermions. 
• Gluonia couplings to 77 

These couplings can be derived by identifying the Euler- 
Heisenberg effective Lagrangian for gg — > 77 via a 
quark constituent loop to the scalar one: £s 77 = 
gSj-ySF$ Fffi , which leads to the sum rule ["]: 

9Sl7 ^ ^V2f s Ml J2 Q\IK^ ( 34 ) 

q—u,d,s 

where Q q is the quark charge in units of e; M Uy d ~ M p /2 
and M s w M^/2; are constituent masses; S refers to 
gluonium (<tb, ...). Then, one predicts the couplings: 

5<t B 77 ~ 9<t' b ti - 5G77 - (°- 4 ~ °- 7 ) a GeV^ 1 , (35) 
and the corresponding widths El : 

12 



o- B ^77 



\9<?Bl~f\ 

16tt 



■M% ~ (0.2 - 0.6) keV 



(36) 



For a self-consistency check of the previous results, one 
can introduce their values into the sum rule: 



(37) 



4 E 5577^2/5 = |^ , 

S'=cr B ,... 

where i? = 3 2Q 2 ., derived [28I3T] from (0|^|7i 72 ), by 
matching the k 2 dependence of its two sides. 

• QCD tests of the <r//o(60Q)= a gluonium ? 
One can notice that the QSSR and LET predictions for 
the mass, hadronic and electromagnetic widths of a low 
mass gluonium as are in remarkable agreement with 
previous results from tttt and 77 scatterings for an on- 
shell or/and Breit-Wigner resonance. However, the fit- 
ted value of direct width in Eq. (24] to 77 needs to be 
improved for a sharper comparison. This "overall agree- 
ment" favours a large gluon component in the tr//o(600) 
wave function. Its large width into tttt indicates a strong 
violation of the OZI rule in this channel and signals large 
non-perturbative effects in its treatment. This large 
hadronic width also disfavours its qq interpretation. In 
fact QSSR predicts for a S2 = {uu + dd), with mass 
of 1 G eV, a tt+tt" width of about (120 - 180) MeV 
[BWpl . and a S 2 -> 77 width of about 5 keV. The 
later can be deduced by using the quark charge squared 



10 However, the analysis in Ref. [4] also indicates that ag having 
a mass below 750 MeV cannot be wide (< 200 MeV) (see also 
some of Ref. 1231 ) due to the sensitivity of the coupling to M a . 
Wide low-mass gluonium has also been obtained using QSSR (1st 
ref. in 1231 ). an effective Lagrangian | 31 | and LET |27l . 

11 This sum rule has been used in 127) for the charm quark. 

12 Due to their M 3 -dependence, the widths of the a' B and G can 
be much larger: (0.4 ~ 2) keV. These widths induce a tiny effect 
of (3-9) xlO" 11 to the muon (g — 2) |66| and cannot be excluded. 

13 The 1st result in 1671 contains an unfortunate numerical error. 



relation (25/9) with r ao _ +77 w 2 keV, obtained from 
the quark triangle loop including the ('(pip) condensate 
contribution [13] . In the same way, QSSR also predicts, 
for a four-quark state, having the same mass of 1 GeV 
|12|13j . a 77 width of about 0.4 eV L 13J. Both qq and 
4-quark scenarios are disfavoured by the value of the 
direct coupling fitted from 77 scattering. 
6. Summary and conclusions 

We have reanalyzed the ttt: and 77 scattering data in 
the mass region below 700 MeV. We applied an analytic 
if -matrix model which is based on an effective theory 
with couplings between resonances, hadrons and pho- 
tons. The model around the threshold region resem- 
bles in its structure with Chiral Perturbation Theory, 
while the inclusion of resonances and a resummation of 
a class of meson loop corrections permit its applications 
towards higher energies. Then, the two-photon decay of 
a resonance can proceed cither through intermediate de- 
cay into hadrons and subsequent annihilation or through 
the direct transition. The low energy processes consid- 
ered here are dominated in our fits by the / (600)/cr 
resonance. Its total radiative width of about 4 keV 
is dominated by the rescattering component which is 
uniquely predicted for known elastic irir scattering. The 
direct coupling width is found an order of magnitude 
smaller (0.13 keV). It can be interpreted as representing 
its intrinsic quark (qq, Aq, . . .) or/and gluon ( gg, . . .) 
constituent contribution to the 77 width. 
A confrontation of the overall fitted resonance parame- 
ters with the QSSR predictions reviewed in Section [5l 
for an unmixed bare gluonium ob having a mass about 
lGeV, indicate that the ct//o(600) can contain a large 
gluon component in its wave function. A light glue- 
ball/gluonium near 1 GeV is also suggested from lat- 
tice QCD calculations including dynamical fermions and 
from the phenomenological analysis of the scalar meson 
spectrum. 

It will be interesting to extend this analysis towards 
higher energies where new channels (KK) will be 
opened and some other resonances will show up. This 
allows for additional constraints on our approach and 
possible tests of the gluon/quark structure of cr//o(600). 
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